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FOREWORD  \ 

'  \ 

The  analysis  presented  in  this  report  was  conducted  when  the  author 
was  connected  with  this  activity.  Since  that  time  the  group  now  engaged 
in  work  on  spherical  shells  have  revised  the  original  manuscript,  par- 
ticularly  the  comparison  of  results  based  on  the  analysis  in  this  report 
with  those  based  on  earlier  analyses. 
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NOTATION 


x,  y,  z 

c< 


ex*  £y 


£i 
£P  J 

e2* 


■ 

t 

O' 

e  , 

v 

Pe  ; 

Ps‘ 

H 

V 

dx»  dy 


♦ 


Coordinates  (see  Figure  2) 
h2/12R2 


Shear  strain 

Strains  in  x-  and  y-directions,  respectively 
EfcCStit^  Strain 
Strain  intensity 
Plastic  strain 

Variational  membrane  strains  in  x-  and  y-directions, 
respectively 

Plasticity  reduction  factor 

e '  » 


Bijlaar  *lasticity  reduction  factor 


> 


reduction  factor 


Gerarc 

Plasticity  reduction  factor  developed  in  present  paper 
Angle  (see  Figure  2) 

Elastic  Poisson1 8  ratio 


I 


Secant  Poisson’s  ratio 
Tangent  Poisson’s  ratio  \ 

Stress  intensity 

Stresses  in  x-  and  y-directions,  respectively 
Shear 

PH  (1  -P82) 

2E^h 


v 


* 


Variational  curvatures  in  x-  and  y- directions, 

respectively 

Angle  (see  Figure  2) 


Primes  refer  to  variational  values  during  the  buckling  process. 


ABSTRACT 


A  solution  for  the  plastic  axi symmetric  buckling  of  thin- 
walled  spheres  under  hydrostatic  pressure  is  derived.  The 
theory  accounts  for  strain -hardening  of  material  and  changes 
of  Poissons  ratio  in  the  plastic  range.  The  plasticity 
reduction  factor  is  expressed  in  terms  of  tangent  and  secant 
moduli  and  new  concepts  of  tangent  and  secant  Poissons 
ratios.  For  typical  engineering  materials  there  is  little 
difference  between  the  results  obtained  from  this  solution 
and  the  earlier  ones  obtained  from  the  solutions  of  Bijlaard 
and  Gerard. 

INTRODUCTION 

Spherical  shells  have  become  more  prominent  in  the  development  of 
submarines  and  deep-sea  vehicles.  Complete  spheres  are  being  used 
for  oceanographic  research  vehicles^  the  TRIESTE  is  a  notable  example. 
On  the  other  hand,  hemispheres  are  being  used  to  effect  closure  at  bow 
and  stern  of  advanced  submarine  designs.  Indications  are  that  spherical 
shells  will  be  used  more  and  more  extensively  in  the  future. 

I 

A  number  of  investigators  have  treated  the  small-deflection  analysis 
of  spherical  shells  subjected  to  external  pressure.  Both  elastic  and 
inelastic  buckling  have  been  studied.  Timoshenko**  summarized  the 
classical,  linear,  small-deflection  theory  for  the  elastic  buckling 

\  \ 

*  References  are  listed  on  page  17. 


pressure  of  complete  spherical  shells  which  was  first  developed  by 

Zoelly  in  1915.  An  expression  for  the  plastic, small -deflection  buckling 

2  3 

of  spherical  shells  was  first  derived  by  Bijlaard.  Gerard  obtained  an 

identical  expression  using  deformation  theory  of  plasticity  in  which  the 
4 

work  of  Stowell  for  flat  plates  was  extended  to  cylindrical  and  spherical 

shells.  Both  Bijlaard  and  Gerard  in  their  analyses  assumed  Poisson’s 

I 

ratio  to  be  equal  to  a  constant,  1/2.  Gerard,  however,  intuitively  modi- 

,  \ 

fied  his  expression  to  include  a  variable  Poisson’s  ratio. 

In  this  report  theory  is  presented  for  the  plastic,  small-deflection 

buckling  of  a  complete  spherical  shell  of  strain-hardening  material  under 

2  3  \  \ 

external  hydrostatic  pressure.  Whereas  Bijlaard  and  Gerard  assumed 
that  Poisson’s  ratio  in  the  plastic  range  is  n.  constant,  1/2,  this  analysis 
considers  a  variable  Poisson's  ratio.  The  results  of  this  more  rigorous 
analysis  are  compared  with  the  existing  analysis  for  spherical  shells,  A 
material  with  an  assumed  stress- strain  relationship  is  used. 

PLASTIC  BUCKLING  THEORY 
VARIATIONS  OF  FORCES  AND  MOMENTS 

Before  the  equations  of  equilibrium  of  an  element  within  a  spherical 
shell  can  be  established,  the  variations  of  the  forces  and  moments  result¬ 
ing  from  the  buckling  process  must  be  obtained.  The  variational  forces 
and  moments  will  be  obtained  by  extending  the  theory  of  Reference  5  to 
spherical  shells. 


2 


In  brief.  Reference  5  uses  the  deformation  theory  of  plasticity  gener¬ 
alized  for  a  variable  Poisson’s  ratio.  The  plasticity  theory  assumes  a 
monotonically  increasing  stress- strain  curve  of  the  strain-hardening 
type;  see  Figure  1.  It  is  assumed  that  Poisson’s  ratio  is  a  function  of 
the  state  of  stress  and  varies  from  its  elastic  value  pe  ■  1/3  to  an  upper 
limit  of  ji8  «  1/2  for  an  isotropic,  plastically  incompressible  solid;  see 
Figure  1. 


Figure  1  -  Variations  of  Strain  and  Poisson’s  Ratio  with  Stress 
The  biaxial  state  of  stress  and  strain  in  a  shell  is  related  to  the  uni¬ 
axial  state  of  stress  and  strain  for  a  simple  compression  specimen  by 

( 

expressions  for  stress  and  strain  intensities.  For  principal  stresses 
and  strains  these  expressions  are  as  follows: 


3 


Stress  Intensity: 


Strain  Intensity: 


The  expression  for  the  plastic  Poisson's  ratio  is  one  derived  by 
Gerard  and  Wildhorn:^ 

1,1  ,  Es 

»snI  “(2"  V  T 

In  Reference  5  these  plastic  stress-strain  relations  were  used  to 
determine  variational  force s  and  moments  for  the  axlsymmetric  buckling 
of  cylinders.  In  this  presentation,  the  axlsymmetric  buckling  of  spheres 
will  be  analyzed.  The  expression  for  the  variation  in  strain  intensity  of 


4 


I 


Reference  5  can  be  readily  modified  to  apply  to  a  sphere.  This  expression 
with  variations  denoted  by  primes,  was  found  to  be: 


£'  *r 


- - }  f<2  -  p.  )  d  +{2}!  -Do  le* 

-.,2ja  d-H  "  s  x  8  yJ  x 


2(  1  -  u  )A  a. 
ps  o  i 


+  [(2(ls-l)Ox+(2-tls)a3r]^rJ 


[6] 


whe  re 

< 

w* 

i 

W 

1 

< 

1 

[(2-^)  -|l-2^)k]2  -3<l-^> 

[7] 

“IT  1 

°  4(1- jAk2 

s 

d 

For  a  spherical  shell  deforming  uniformly  d  =  d^j  therefore  k  a  »  1 

and  K2  *  1.  Also, 

X 

x” 

II 

i 

N 

X 

i— •- 

[8a] 

ey=  ^-zxi 

1 - 1 

JO 

CO 
i — I 

Equations  [6]  and  [7]  become 


,  ti+e2  (xj  +  x'J 

h  B  2(1  -  \ij "  2(1  -  n  )A 


C9] 


s  o 


s  o 


2(  1  -  u  ) 


Let 

Then 


S  ■  E  e  . 
x  s  x 


S  -  E  e»  -  —  (E  -  E  )e* 
x  s  x  8  t  i 


fio] 

[11] 

[12] 
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Let  z  =  z  when  e!  «=  0;  then,  from  Equation  [9], 

O  1 


£1  +  ^ 


Z°  '  xi  +  x'z 


Substituting  Equations  [8a],  [9],  and  [l 3]  into  Equation  [12]  gives 


sx"  E.(  *1  -  *Xi>  +  2(1*-  u)A  1  (Xi  +  Xi>'*  -  *„>  [H] 


S  O  1 


Similarly, 


e  (E  -  E  ) 


VE.(3-»xi»+gi^ r 

O  X 


<Xi  +  Xi«“  -  =„>  [15] 


S  +(i  S 

Since  o  e  * *>  1 

x  .  2 

1  -  u 
rs 


X  •  s 

0  - - - 

x  ,  2 

1  -  K 


The  variational  axial  force  is 


SX  +  -SS;  MSV  ■'2*13Sx]l1S 


(1  - It 


^  r+1 

N*°J 


d'  dz 
X 


Substituting  Equation  [16]  into  Equation  [l?]  and  integrating  gives 


I,x-Bp,AlH  +  '*.Al2^ 


[18] 


where 


6 


,  (l-VsVw 

"  *  “  <  >  f  B  I 


4A  ( 1  -  u  ) 
o 


j  i1  -  Vs)'1  +  ■>,». 


4u  A  a  -  U  ) 


[21] 


Similarly, 


The  variational  axial  moment  is 


f+h/2 

u  / 


[23] 


Substituting  Equation'll]  into  Equation  [23]  and  integrating, 

where 


D  = 


E  h 
s 


P  1- 


[25] 


Similarly, 


The  plasticity  coefficients  A^  and  reduce  to  expressions  obtained 


by  Gerard  when  is  set  equal  to  1/2. 
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EQUATIONS  OF  EQUILIBRIUM 

The  equations  of  equilibrium  for  a  shell  buckling  axisymmetrically  are 
given  by  Timoshenko*  as:* 


dN 


t  r 


»  PR 


+(N  -N  )  cotfl  -  Q  -  , 
d0  x  y  x  2 


(u  dw  \ 
R  +  Rd9  / 


dQ 


i  i 


( 


du 


-jr-  +  Q  cotO  +N  +N  +  PR  I  ~  +~cotO 
d0  x  x  y  \  Rde  R 


PR  /  du  d2w  \  PR  _  /  u  dw  \ 

'Tl555+^/'~co8U+535) 


dM  ,  ,  , 

+  (M  -  M  )  cotO  -  Q  R  *  0 
d6  x  y  x 


2w\ 
“  »  / 


[27] 


The  notation  used  is  explained  in  Figure  2.  -n  terms  of  the  displacements,* 
the  variational  membrane  strains  and  variational  curvatures  are:* 


,  du  w 
el"  RdO  "  R 


,  u  w 

£2  *  R  Cot  ^  ”  R 


Xl“ 


.2 

d  w 


du 


R2d02  R2d0 


i  / u  dw  \  0  v 


[28] 


can  be  eliminated  by  combining  the  first  two  qualities  of  Equations  [27] 
By  substituting  Equations  [18],  [22],  [24],  [26],  and  [28]  into  the 
resulting  two  equations,  there  results 


*  The  displacements  u  and  w  are  variational  values  but  for  convenience 
will  not  b<»  denoted  by  primes. 


I 


u  +  0<,[Ax.~f +  Aicot855--<i‘1 


sA12  +  A1  c°‘  81  “ 


"  (A1  +  l*SA 


12 1  df  +  «  [A1  +  A1  cot  8  *  ( ^aA12  +  A1  c»t2  8  >  if  ] 


'  *  (“  +df  )  ' 


li  +i*.Ai2,(; 


+  ucotO  -  2w 


)  +  «[-  A  -[■-  2A  cotQ  ^4 

’  L  deJ  de2 


o  « 

+  (A1  +  ^12  +  Alcot  0)  dO  "  COt0  (2A1  "  ^12  +A1cot20)  u 


A1  -  2AX  cote  +(AX  +  n  A  +  A  cot20)  ^4 

1  de  1  _  de3  1  8  12  1  d02 


•cot8<2A1-  %A12  +Alcot29)|:j 


wher* 


<t>  I  -  a  cot  e  -  ^  +  4w  +  cot  e 


°<  *  h^l2B2 


dw  ,  d2w  \ 

dT+^J 


PB(1  -  vT) 
*  - - L£_ 

2E  h 


I 


Equations  [29]  and  [30]  reduce  to  Equations  [e]  and  [f]on  page  492  of 
Reference  1  when  L'  =  E  and  p  =  p  . 


Figure  2  -  Forces  and  Moments  on  an  Element  of  a  Spherical  Shell 

PLATE  BUCKLING  EQUATION  I 

\ 

\ 

The  solution  for  the  buckling  pressure  is  practically  identical  to  that 

i 

described  by  Timoshenko  in  Reference  1.  Thus,  introducing  a  new 
variable,  the  angle  ^  (see  Figure  2)  as  follows: 

a  = 

d9 

and  using  the  symbol  H  for  the  operation 
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there  is  obtained  from  Equation  [29] (after  integrating  once)  and  also 
from  Equation  [30]  the  following:  I 

\ 

AjHCiji)  +  <AjH{w)  -  (Aj^  +  ^sa12)(^  +  w)  -  <fr  {$  +  w)  =  0  [33] 

o(A ^HH( +  w)  -  (A1  +  \ifiAlz)  H(^>)  -  (3AX  +  ji^A^)  «H(w  +  $) 

+  2(A^  +  ji  A^2)(t|i  +  w)  +  $  [-H(iJ>)  +  H(w)  +  2($  +  w)]  »0  ''  [34] 

As  in  the  case  of  the  elastic  theory.  Equations  [33]  and  [34]  have 
solutions  in  terms  of  spherical  functions.  Bather  than  repeat  steps 
previously  described,  *  only  the  final  buckling  equation,  resulting  from 
setting  the  determinant  of  two  homogeneous  equations  equal  to  zero,  is 
‘  presented: 

-  W 

EquationX3S^*reduces  to  the  classical,  elastic,  small-deflection 
theory  Eg  *  E  and  =  jig. 

The  plastic  buckling  equation  is  more  easily  related  to  the  elastic 
'  solution  if  the  concepts  of  secant  and  tangent  Poisson's  ratios  are  intro- 
■  disced.  The  secant  Poisson's  ratio  is  the  weighted  average  of  the  elastic 
Poisson's  ratio  ^  and  the  fully  plastic  value  see  Figure  1.  Thus 


But  £e  =  C'J  E  and  £  =  c^E  ,  thus 


s  s 


H 


Since  e  s  e  -  £  ,  then 
p  i  e 


p  ’  %  -  ee  "  '  _1_  I 

Es  E  E»  E 


J_  _1_  1 

E  *  E  "  E 
P  « 


Substituting  Equation  [38]  into  Equation  [37], 


P  =  p  -  {  U  -  11  )  TT 
s  rp  rp  re  E 


For  a  fully  plastic,  isotropic,  incompressible  material,  p^  =  1/2.  Then, 

E 

-  1  Z1  1  S 

vr'rvr  !  [40] 

i 

Equation  [40]  is  identical  to  Equation  [5]  presented  by  Gerard  and 
Wildhorn^  as  the  plastic  Poisson’s  ratio. 

The  tangent  Poisson’s  ratio  is  defined  herein  as  the  number  which, 
when  multiplied  by  the  variation  in  strain,  gives  the  variation  in  strain 
In  the  transverse  direction.  Thus 


l*tei  "  ’•‘.V  =  Eil*i  +  S  *1 


I 


-Therefore, 


' - 4 - 


The  variation  of  Equation  [40]  is  I 

\ 

1  *i 

*‘«,^(Es/E-Et/E)  t 

Substituting  Equations  [40]  and  [43]  into  Equation  [42]  gives  the  final 


expression: 


“t  2  (?-  ‘‘e’r  - 

In  terms  ofTthe  secant  and  tangent  Poisson’s  ratios.  Equations  [39] 
and  [44],  Equation  [35]  can  be  show  to  be: 


VE.  M 


P  =2 
cr 


DISCUSSION  AND  CONCLUSIONS 

Equation  [45]  may  be  compared  with  the  classical,  elastic  solution 


P  = 
e 


'  3(  1  -  u  ) 


through  the  use  of  the  plasticity-reduction  factor  rj  defined  as 


i 


Substituting  Equations  [45]  and  [46]  into  Equation  [47]yields 


E  ^  Es(l  -  ,it)(l  +  jis)  L” 

Equation  [48]  reduces  identically  to  that  developed  by  Bijlaard2  and 

3 

Gerard  if  p  and  p  are  set  equal  to  1/2;  thus 

t  -  s _ _ 


E.ll  -  ,  1 


E 

.  - 
^B,G  ”  E  1 


Et(i  - »: 

E  (0.75) 
s 


3  . 

Equation  [48]  also  reduces  to  that  proposed  by  Gerard  when  he  intui- 

♦ 

tively  accounted  for  a  variable  Poisson's  ratio  in  the  plastic  range,  if  p^ 
is  set  equal  top^;  thus 


nG  *  E  „ 


r.d  -  p6) 
t  re 

E  (1  -  p2) 
s  r  s 


Equation  [45],  therefore,  reduces  to  the  following: 

1.  Zoelly's  classical,  small  deflection,  elastic  theory  when  moduli 
and  plastic  Poisson's  ratios  take  on  elastic  values;  see  Equation  [48], 

2.  Bijlaard's  and  Gerard's  plastic  solution  when  secant  and  tangent 

! 

Poisson's  ratios  are  set  equal  to  1/2;  see  Equation  [49]. 

3.  Gerard's  intuitively  modified  plastic  solution  when  the  tangent 
Poisson's  ratio  is  set  equal  to  the  secant  Poisson's  ratio;  see  Equation 
[50]. 
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The  collapse  pressures  as  determined  by  the  inelastic  theories  of 
Bijlaard  and  Gerard  are  compared  with  Equation  [45]  of  this  paper  in 
Figure  3.  The  abscissa  is  the  ratio  of  the  classical,  small-deflection, 
elastic  collapse  pressure  to  the  pressure  at  which  the  average  stress 
reaches  the  yield  point,  P^,  as  determined  by  the  0,2  percent  offset 
method.  The  ordinate  is  the  ratio  of  theoretical  inelastic  collapse 
pressure  Pcr,  (as  determined  by  Equation  [45]  ,  Bijlaard  and  Gerard)f to 
Py,.  A  typical  stress-strain  curve  for  7075-T6  aluminum  bar  stock  has 
been  assumed  in  all  theoretical  calculations.  There  is  little  difference 
between  the  pressures  obtained  from  the  inelastic  theories  of  Bijlaard 
and  Gerard  and  the  more  rigorous  theory  developed  herein.  The  maxi¬ 
mum  difference  between  the  three  theoretical  collapse  pressures  for  a 
spherical  shell  of  7075-T6  aluminum  is  about  2  percent;  see  Figure  3. 
Similar  results  would  be  obtained  for  other  strain-hardening  materials. 

Theory  which  fails  to  consider  the  effects  of  imperfections,  residual 
stresses,  boundary  conditions,  and  penetrations  will  not  consistently 
predict  the  collapse  strength  of  spherical  shells  normally  encountered 
in  engineering  practice.  However,  this  theory  will,  if  verified  by 
experiment,  serve  as  a  reference  to  which  tests  of  fabricated,  spherical 
shells  may-bfc’TTaTripared  to  determine  the  detrimental  effects  of  initial 
departures  from  sphericity,  variations  in  shell  thickness,  residual 
stresses,  boundary  conditions,  and  penetrations. 
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